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Abstract 



By extending the Standard Model with three right-handed neutrinos 
{Ni) and a second Higgs doublet {H2), odd under a Z2 symmetry, it is 
possible to explain non-zero neutrino masses and to account for the dark 

CD ■ 

matter. We consider the case where the dark matter is a scalar and study 
its coannihilations with the right-handed neutrinos. These coannihilations 
tend to increase, rather than reduce, the dark matter density and they 
^N) I modify in a significant way the viable parameter space of the model. In 

^ ■ particular, they allow to satisfy the relic density constraint for dark matter 

masses well below 500 GeV. The dependence of the relic density on the 
relevant parameters of the model, such as the dark matter mass, the mass 
splitting, and the number of coannihilating fermions, is analyzed in detail. 
We also investigate, via a scan over the parameter space, the new viable 
^NJ , regions that are obtained when coannihilations are taken into account. 

' Notably, they feature large indirect detection rates, with {gv) reaching 

CO ' values of order 10~^''cm^s~^. Finally, we emphasize that coannihilation 

effects analogous to those discussed here can be used to reconcile a thermal 
freeze-out with a large {av) also in other models of dark matter. 

•rH ' 

X 

H ■ 1 Motivation 

■ 

Dark matter and neutrino masses provide the only two clear evidences we cur- 
rently have of physics beyond the Standard Model (SM). Cosmological observa- 
tions indicate that about 25% of the energy density of the Universe consists of 
a new form of matter usually called dark matter [T]. Since none of the known 
elementary particles satisfy the conditions required to explain the dark matter, 
new physics should be responsible for it. Even if we still have no clue of what 
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this new physics is, circumstantial evidence suggests that it lies at the TeV scale 
[5] . The observation of neutrino oscillations [31 HI [5J |5] implies that neutrinos 
have non-zero masses -a fact that cannot be accounted for within the SM. In 
extensions of the SM, neutrino masses can be generated in several ways. Some of 
them require new physics at very large scales whereas others do not. Normally, 
dark matter and neutrino masses are treated as separate problems, without any 
relation to each other. It may be, though, that they are intimately linked and 
that both originate from physics at the electroweak scale. 

A simple model that explicitly realizes this idea was proposed several years 
ago in [7 . In this model, known as the radiative seesaw model, the SM is 
extended with a second Higgs doublet and three right-handed neutrinos. A 
viable dark matter candidate is obtained by assuming that the new fields are 
odd under a Zi discrete symmetry, and neutrino masses are generated by one- 
loop processes mediated by the odd fields, in particular, by the dark matter 
candidate. The dark matter phenomenology of this model has been extensively 
studied in previous works, see e.g. [H HI HOl El [H US HI [HI [H] . None of these, 
however, has considered the possible effect of coannihilations with the right- 
handed neutrinos for scalar dark matter. That is precisely what we do in this 
paper. We will see that these coannihilations tend to increase the relic density of 
dark matter and may modify in a significant way the viable parameter space of 
the model. Specifically, they allow to reduce the lower bound on the mass of the 
dark matter particle from 500 GeV down to about 100 GeV. We study in detail 
how these coannihilation processes affect the calculation of the relic density for 
different values of the parameters of the model, and we analyze the resulting 
viable parameter space. Remarkably, we find that, due to coannihilations, the 
indirect detection rate of dark matter is strongly enhanced within the new viable 
regions. 

The rest of the paper is organized as follows. The next section introduces 
the radiative seesaw model and fixes our notation. In section [31 we discuss 
the coannihilation formalism and how it can be applied to our model. Our 
first results are presented in section 01 where we analyze the dependence of the 
coannihilation effects on the dark matter density for different sets of parameters. 
In section[5lwe investigate the new viable parameter space that is obtained when 
coannihilations are taken into account. In particular, we demonstrate that it 
features much larger values of (crw) . Section[Slbrings up some possible extensions 
of our work and section [71 presents our conclusions. 

2 The radiative seesaw model 

The radiative seesaw model is a minimal extension of the Standard Model 
that simultaneously accounts for neutrino masses and for the dark matter. It 
includes three right-handed neutrinos Ni {i = 1,2,3) and one additional Higgs 
doublet H2- These new fields are assumed to be odd under a new (and exact) 
Z2 discrete symmetry under which all the SM fields are even. The new terms 
that appear in the Lagrangian of this model are 

Cn = N,i0PRN, + (D^H^)^ {D^'H2) - ^N^PrN, -f h^J^HlPRN, + h.c, (1) 
and 
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-Lh.m. = lAaliii + M2^]^2 + y + Y (^^1^2 j (2) 

+ A3 (b\b^ + A4 (i/Ii/2) (i^l^fi) + Y (ir|^f2)' + h.c, 

where il\ is the SM Higgs doublet. Had only H2 (but not Ni) been added 
we would have obtained the so-called inert doublet model [T71 [TB] , a minimal 
extension of the Standard Model that can explain the dark matter. Thus, in 
certain regions of the parameter space, the radiative seesaw can be seen as an 
extension of the inert doublet model that accounts also for neutrino masses. 
It is also important to notice that A^i, being odd under Z2, are not the same 
right-handed neutrinos that appear in the usual seesaw mechanism. They do 
not give rise, for instance, to a Dirac mass term for the neutrinos. 

The new scalar sector of this model contains four physical states: two charged 
states, iJ^, and two neutral ones, and Their masses are given by 

M^o = ^2 _^ (A3 + A4 + A5)^;^ 

Mlo -A^^ + (A3+A4-A5)^;^ (3) 

where v is the vacuum expectation value of iJi . We take as the free parameters 
of the scalar sector the physical masses (Af//o,M^o,M^±) and the coupling 
A = A3 -|- A4 -|- A5 . 

Majorana neutrino masses in this model are generated at one-loop by the 
exchange of ° , and Ni . The resulting mass matrix is given by 



where Mi are the masses of the right-handed neutrinos, Mq ~ /ij + (A3 + X4)v'^, 
and the loop function is given by 

\ 2; / xlogx\ 
I{x) = 1 + . (5) 



1 — a; V 1 — X 



Compatibility with present neutrino data can be achieved if the matrix of 
Yukawa couplings hai has the following structure [TT] 



(6) 



From it, it follows that 6*23 = 7r/4, 6*13 = and tan = -^h'^/h^. Current 
data, see e.g. [1^], then implies h'^/h^ « 0.95. The relation between neutrino 
masses and the parameters hi is as follows 




m: 
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where 



) 



(8) 



and we have taken hi to be real. Notice, therefore, that only one parameter, 
from those appearing in equation ([6]) , remains undetermined by the experimental 
data. We take that free parameter to be /i2. hi and /13 are then calculated, 
for given masses of the odd particles, from equations (O and ([5]). In any case, 
as we will see later, our results do not depend on the specific choice made in 
equation 

In this model, one-loop diagrams mediated by the odd particles give rise 
to lepton flavor violating processes, such as /i — 67 and r fj.'j, which may 
impose important restrictions on the parameter space of the model. It turns 
out, though, that in the region of interest to us, the Yukawa couplings are small 
and the experimental bounds on lepton flavor violation are easily satisfied. 

Regarding dark matter, in this model the lightest odd particle can either 
be a fermion (one of the three right-handed neutrinos) or a scalar or A"). 
The former case has been studied in [SlinilTnilllllllllHIIl] and the latter in 
[TSIITB]. In this paper, we focus on the case where the dark matter candidate is a 
scalar, H'^ . Naively, one might expect the resulting dark matter phenomenology 
to be similar to that of the inert doublet model, but that is not necessarily 
true. As we will show, H^-Ni coannihilations may play an important role in 
the determination of the relic density, modifying the viable regions and the 
prospects for the detection of dark matter. 

3 Coannihilations and their effect on the rehc 
density 

The observed dark matter density is believed to be the result of a freeze-out 
in the early Universe. This freeze-out process can be modified by the presence 
of additional particles that are close in mass to the dark matter particle -what 
is usually known as coannihilations [20j . Coannihilations may lead both to an 
increase or a decrease of the final abundance of the lightest stable particle, the 
dark matter candidate. In the context of supersymmetry, for example, coannihi- 
lations usually cause a suppression of the relic abundance. Since the ncutralino 
is bino-like in most of these models, it annihilates rather inefficiently in the 
early Universe and for that reason is typically too high. It turns out though 
that, in certain regions of the parameter space, the ncutralino becomes almost 
degenerate with a stau so that stau-neutralino coannihilations become relevant. 
Because staus annihilate more efficiently than bino-like neutralinos, the effec- 
tive coannihilation cross section is larger than the ncutralino annihilation cross 
section and the final relic density is lower than without stau coannihilations, 
making it easier to satisfy the dark matter bound. 

The opposite behavior is observed in scenarios with Universal Extra Di- 
mensions (UED) ,21j. In UED, the spectrum of Kaluza-Klcin states is highly 
degenerate so coannihilations are expected to play an important role. It has 
been shown that, for realistic spectra, the ultimate effect of these coannihila- 
tions is to increase the relic abundance of the lightest Kaluza-Klein particle, 
reducing the dark matter mass required to satisfy the relic density constraint 
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[321 [231 [23 • Within supersymmetry, it is also possible to find examples where 
the effect of coannihilations is to enhance rather than reduce the relic density. In 
[25| , it was shown that if bino-like neutralinos annihilate through the resonant 
exchange of a heavy Higgs, slepton coannihilations give rise to an increase in the 
relic density. In [26], models where the neutralino is higgsino-like or wino-like 
were considered and it was shown that slepton coannihilations not only lead 
to an increase in the relic density but also to an enhancement in the predicted 
indirect detection signals. 

Let us now briefly review, following |201 126] , the coannihilation formalism 
and use it to obtain an estimate of the expected effect within the radiative seesaw 
model. The effective annihilation cross section for a system of n coannihilating 
particles is given by 

n 

where — {mi — m^)/my^ is the relative mass splitting between the «-th coan- 
nihilating particle and the lightest stable particle X; a; = m^/T, Uij denotes the 
annihilation cross section of particles i and j into Standard Model particles, gi 
stands for the number of internal degrees of freedom of particle i and 

n 

ffeff = ^5,:(l + A,)3/2e--^'. (10) 

1=1 

Here, we are interested in the case where the dark matter particle, , has an ef- 
fective pair annihilation cross section much larger than that of its coannihilating 
partners, the right-handed neutrinos Ni. This situation arises naturally within 
the radiative seesaw model. In fact, annihilates efficiently into W'^W" and 
Z'^ Z'^ via gauge interactions whereas the interactions of the right-handed neu- 
trinos are determined by the Yukawa couplings and can easily be suppressed. It 
must be kept in mind, however, that the H'^ effective annihilation cross section 
without Ni coannihilations is actually the result of coannihilations within the 
scalar sector. As we will see, the mass splitting of and with H'^ is al- 
ways small within the viable regions, so those coannihilations are unavoidable. 
If we denote by (Jh^h" = '^cS effective H'^ annihilation cross section with- 
out right-handed neutrino coannihilations then we have that aj^ojjo ^ (^H°Ni^ 
<^NiNi, where <JH°Ni^ (^NiNi respectively indicate the Ni coannihilation and self- 
annihilation effective cross sections. Denoting with the effective degrees of 
freedom when the Ni are much heavier than , the new effective total annihi- 
lation cross section cr^^ can be written as a function of the effective number of 
degrees of freedom g^L including the Ni as 



^Ni „ / 5cff(,^f.o. 



<ff - ^cff ^^f™ (11) 

\9^s[xi.o.) J 

where Xf.o. corresponds to temperatures around the freeze-out. Since g^ff{xf,o.) 
is larger than (jiJ?g(a;f.o.), the effective annihilation cross section decreases with 
the inclusion of Ni coannihilations -and the dark matter density increases ac- 
cordingly. In the limit where M^o — MjV; becomes negligible (A^ — 0) one can 
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Mj^o [GeV] 



Figure 1: The relic density as a function of Mho for different values of Mn. — 
Mffo. In this figure we have set A = 0.01, /12 = 0.01, M^o — Mu± — M[j« + 
5 GeV and we have assumed that the three ferniions have the same mass: M^-^ — 
MjVj = MjVg . Notice that fi/i^ decreases with increasing M^. — Mho . 



write 



(12) 



for the ratio between the relic density including Ni coannihilations (fi^*) and 
that without doing so {^^). This approximate expression depends only on 
the internal degrees of freedom that take part in the coannihilation process. 
In our case, there are four scalar degrees of freedom {H^, A^, H^) so go — 
4. Each right-handed neutrino has two degrees of freedom and we could have 
coannihilations with 1, 2 or 3 of them, so gN- can respectively take the values 2, 4 
or 6. The ratio fi^'/fi'^ is then equal to 9/4 for 1 coannihilating neutrino, 4 for 
two coannihilating neutrinos, and 25/4 for 3 coannihilating neutrinos. These 
numbers, being obtained for a negligible mass splitting between Ni and H'^, 
actually represent an upper bound to this ratio, so we expect to find smaller 
effects in more realistic situations. From this simple discussion we can already 
conclude, in any case, that coannihilations with right-handed neutrinos can 
increase the relic density by up to a factor 6 or so. 

In the next section we evaluate numerically the effect of coannihilations on 
the relic density. To that end, we have implemented the radiative seesaw model 
into micrOMEGAs [37], which automatically takes into account all possible 
coannihilation effects. 
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Figure 2: The relic density as a function of M/^o for different values of 
MmJMho- In this figure we have set A = 0.01, /i2 = 0.01, M^o = Mu± = 
Mfjo +5 GeV and we have assumed that the three fermions have the same mass: 
Mat^ = Matj = Matj. Notice that U.h'^ decreases with increasing MjqjMjjo. 



4 Results for the relic density 

In this section, the effect of Ni coannihilations on the relic density of is 
numerically studied. Specifically, we analyze how ilh^ depends on the Ni-H° 
mass splitting, on the value of the neutrino Yukawa couplings and on the 
number of coannihilating neutrinos. We also determine some viable regions and 
study their variation with the parameters of the model. 

When coannihilations with A^; are negligible, the dark matter phenomenol- 
ogy of this model reduces to that of the well-known inert doublet model in the 
heavy mass regime [HI [IH]- It is useful, therefore, to keep in mind the main 
two features of this latter model. First, the mass splitting between the scalars 
is necessarily very small along the viable regions. Second, the viable mass range 
starts at Mj^o > 500 GeV. For smaller values of Mfjo, Ufjo is always below the 
observed range. We will demonstrate that these features are modified in the 
presence of coannihilations with the right-handed neutrinos. 

Figured] shows the relic density as a function of Af^o for several values of the 
mass difference: M^i — Mho — 1 GeV (dashed line), 3 GeV (solid line), 10 GeV 
(dash-dotted line). For comparison, the case where coannihilation effects are 
negligible, A/jy. ^ M//o, is also shown (dotted line). The horizontal band 
shows the WMAP range for Q,h^ [T] . We see that coannihilation effects increase 
the relic density and consequently reduce the value of M[jo that is compatible 
with WMAP data. As the mass splitting is reduced, coannihilation effects 
become more relevant and the increase in the relic density is larger. The viable 
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Mjjo [GeV] 

Figure 3: The relic density as a function of M^o for different values of In 
this figure we have set A = 0.01, Mj^a — Mh± — M^o + 5 GeV and we have 
assumed that the three fermions have the same mass: Mj^-^ — Mj^^ — Mj^^ — 
Mho + 5 GeV. 



value of Mffo, for instance, goes from about 750 GeV for M^. ^ Mho (no 
coannihilations) to about 250 GeV for M^^. — Mho = 1 GeV. 

To better understand the degree of degeneracy between Ni and Mho that 
is required to obtain a certain effect on the relic density, it is useful to look at 
the dependence of ^h? with the relative mass difference, [Mm^ — Mho) /Mho 
-see figured] Notice that a 10% degeneracy (dash-dotted line) produces almost 
no visible effect, and that in that case Mho is compatible with the observations 
for Mho ~ 750 GeV -the same value obtained without coannihilations (dotted 
line). A 1% degeneracy (solid line) reduces the viable value of Mho to about 
400 GeV, and a 0.1% degeneracy lowers it further to about 150 GeV. 

The neutrino Yukawa couplings also affect the relic density, as they determine 
the strength of the H^-Ni cross section. We expect, in particular, that only for 
small yukawas the right-handed neutrinos will behave as additional degrees of 
freedom and increase the relic density. Figure [3] shows ^h? as a function of 
Mho for different values of /i2. For /12 = 0.1 (dotted line) we notice that the 
relic density is always below the WMAP bound. For /12 = 0.05 (dash-dotted 
line) Vlh? is slightly larger and it becomes possible to satisfy the dark matter 
constraint in the range 500 < Mho /GeV < 600. For /12 = 0.01 (solid line) the 
viable value of Mho is smaller, about 400 GeV. That same value is obtained 
for /i2 = 0.001 (dashed line). In fact, these last two lines only deviate from 
each other at large values of Mho, where fi/i^ is well above the observed range. 
Considering even smaller values of /12 would have no further effects on the relic 
density. 
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M^. = Mjj. = M^i + 5 GeV 




Figure 4: The relic density as a function of M^o for different numbers of coan- 
nihilating N. In this figure we have set A = 0.01, /i2 = 0.01, M^o = M^± = 
Mho + 5 GeV, and we have assumed that the coannihilating particles have a 
negligible mass splitting: M^. — M^o ^ Mj^o . 



In all the previous figures we have assumed that the three right-handed 
neutrinos have the same mass, M^-^ — M^^ = M^^, so they all coannihilate 
simultaneously with . It may well be, though, that only 1 or 2 of them actu- 
ally take part in the coannihilation process. In the next figures, we study how 
the predicted relic density depends on the number of coannihilating N . Figure 
m shows the relic density as a function of Mh« for 0, 1, 2, and 3 coannihilating 
fermions. In this figure we have assumed that the mass splitting between the 
coannihilating neutrinos and is negligible. Ma?. — M^a <^ M^^o, whereas the 
mass of the non-coannihilating neutrino is much larger than that of H'^ . First 
of all, notice that, as expected, the relic density increases with the number of 
coannihilating particles or rather with the number of coannihilating degrees of 
freedom. If there are coannihilating neutrinos (dashed line) we reproduce the 
known result from the inert doublet model and we observe that M^o should be 
around 750 GeV to account for the dark matter. For 1 coannihilating neutrino, 
the viable mass moves down to M^o ~ 300 GeV, a value that can never be 
reached within the inert doublet model. For 2 and 3 coannihilating neutrinos 
we obtain respectively Mh« ^ 200 GeV and M^o ~ 150 GeV for the values 
that are compatible with WMAP. Thus, the number of coannihilating neutrinos 
plays a very important role in the determination of the i?" relic density. 

Figure [S] shows flh^ as a function of the relative mass difference, {M^. — 
Mfjo)/M}jo, for all the possible values of the coannihilating neutrinos and a 
fixed value of the dark matter mass. Mho = 350 GeV. First of all notice that 
the coannihilation effect completely disappears when (A/at. ~Mho)/Mh« reaches 
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Figure 5: The relic density as a function of the ratio {Mn. — Af^o)/Af^o for 
different numbers of coannihilating N. In this figure we have set A — 0.01, 
/i2 = 0.01, Mao = Mh± = Mho + 1 GeV, and Mho = 350 GeV. 



about 10%. For higher values, the relic density is not affected by coannihilations 
so it becomes independent of the number of coannihilating N. It is also clear 
from this figure that flh^ cannot be made arbitrarily large by reducing the mass 
splitting. There is a limit to this trend, for Qh^ reaches an asymptotic value at 
very small {Mn.— Mho) /Mho. In that limit, the numerical results coincide with 
the estimates obtained in section[3] For coannihilating neutrinos (dashed line) 
we observe from figure O that the dark matter constraint cannot be satisfied -as 
expected, because that is the inert Higgs model limit and Mho — 350 GeV < 
500 GeV. For 1 coannihilating neutrino, the right relic density is obtained for 
a mass degeneracy below the 0.1% level. The required mass degeneracy lies 
between 0.5% and 1% for 2 coannihilating neutrinos, and between 1% and 2% 
for 3 coannihilating neutrinos. 

It is also important to look at how the viable regions, those that satisfy the 
dark matter constraint, depend on the mass splitting and on the number of 
coannihilating particles. Figure [5] shows such regions in the plane {Mho ,{Mf^. — 
Mho)/Mho) for 1, 2 and 3 coannihilating neutrinos. Notice that these regions 
start respectively at about 150 GcV, 200 GeV, and 350 GeV, and that in all cases 
the initial degeneracy is below 0.1%. As Mho increases so does the required 
degree of degeneracy, reaching about 1% respectively for 350 GeV, 400 GeV, and 
550 GeV. When [Mm^ — Mho) /Mho becomes of order 10% the coannihilation 
effect fades away and all three lines converge. 
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Figure 6 : The regions in the plane {M[jo , {Mjy. — Mua ) j M^a ) that are consistent 
with the rehc density constraint for different numbers of coannihilating iV. In 
this figure we have set A = 0.01, hi = 0.01, M^o = M^i = M//o + 5 GeV. 



5 Implications 

We have already observed that H'^-Ni coannihilations can modify the predicted 
relic density of iJ", giving rise to new viable regions which may differ consider- 
ably from those obtained without coannihilations. In this section we are going 
to study such viable regions in more detail. We will see that they indeed allow 
for significant differences in the spectrum and result in modified prospects for 
the detection of dark matter. 

In this analysis, we will make use of a scan over the entire parameter space of 
this model. We have generated four large sample of models, each with a different 
number of coannihilating N. In each case, we have varied the parameters of the 
model within the following ranges: 

100 GeV < Mho < 1 TeV (13) 

Mho < Mao < Mho + 40 GeV (14) 

Mho < Mh± < Mho + 40 GeV (15) 

Mho < Mn, < Mho + 40 GeV (16) 

10"^ < A < 10-^ (17) 

10"^ < /i2 < 10"^ (18) 

Here, M^r. denotes the masses of the coannihilating N. The non-coannihilating 
neutrinos were given a mass much larger than Mho . From the randomly gener- 
ated models, we have selected those that satisfy the dark matter constraint to 
obtain the sample of viable models that we are going to analyze next. 
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Figure 7: A scatter plot showing the mass difference between the scalar particles 
versus Af^o for different numbers of coannihilating N: 3 (red), 1 (blue) and 
(black). 




10 20 30 40 10 20 30 40 



Figure 8: Scatter plots showing the mass difference between the scalar particles 
versus M^. — M^o for two different numbers of coannihilating iV: 3 (red), 1 
(blue). 



Figure [7] shows a scatter plot of the mass difference between the scalar parti- 
cles versus M[jo for different numbers of coannihilating N: 3 (red), 1 (blue) and 
(black). The left panel displays Mao — Mjfo and the right panel M[j± — Mjjo. 
When there is no coannihilations (black points), we notice that there are no 
points below Mj^o ^ 500 GeV and that the mass splitting is always very small, 
< 10 GeV. These are the well-known results for the inert doublet model. If 
there is one coannihilating neutrino (blue circles), the viable mass range starts 
at about M^a ^ 300 GeV and the mass splitting can be larger, reaching 
Myio — Mfjo ^ 25 GeV and Mh± — M^o ^ 20 GeV. In the case where the three 
neutrinos coannihilate (red squares), the whole range of M^o becomes viable 
and the mass splitting can reach the maximum values we examined, ^ 40 GeV. 
It is clear, therefore, that the viable regions when H^-Ni coannihilations are in- 
cluded may feature much smaller dark matter masses (down to A/^o ^ 100 GeV 
or so) and larger values for the mass splitting between the new scalar states. 

What allows Mao— Mho and Af^± — Mho to be larger are the coannihilations 
with the right-handed neutrinos, as illustrated by figure |8l It displays the scalar 
mass differences versus M^. — A/^o for 1 (blue circles) and 3 (red squares) 
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Figure 9: Scatter plot of the spin-independent direct detection cross section (crsi) 
versus Mfjo for (blue circles) and 3 (red squares) coannihilating neutrinos. In 
this figure we show only those models featuring dsi > 10^^^ pb at tree-level. 
The solid line corresponds to the present bound from XENONIOO, while the 
dashed line corresponds to the expected sensitivity of XENONIT. 



coannihilating neutrinos. Notice, from both panels, that Mao —Mho and Mh± ~ 
Mfjo can become relatively large only when Mjy. — M^o becomes small. Thus, 
one way or another there must always be some particle highly degenerate with 
Mho. 

Let us now look at how the prospects for the detection of dark matter are 
modified in the presence of coannihilations. The dark matter spin-independent 
direct detection cross section (crsi) in this model is determined at tree-level by a 
Higgs mediated diagram and is proportional to A^-just as for the inert doublet 
model. As pointed out recently in ^29^, however, the one- loop corrections to this 
cross section can be large and substantially modify the predicted value of asi- In 
particular, they always bring its value within the reach of future direct detection 
experiments, ~ 10~^^ pb. These results are not affected by the presence of 
coannihilating right-handed neutrinos. Their main effect is simply to allow for 
viable models with smaller masses for the dark matter particle, a region where 
direct detection experiments have a larger sensitivity. This fact is illustrated in 
figure [21 which shows a scatter plot of crgi versus Mho for (blue circles) and 3 
(red squares) coannihilating neutrinos. The solid line shows the present bound 
from XENONIOO [30] and the dashed line the expected sensitivity of XENONIT. 
Notice that whereas for 3 coannihilating neutrinos a handful of models with 
masses around 200 GeV is already ruled out by the experimental bound, not 
a single model is currently above the limit for coannihilating neutrinos-they 
all lie at least one order of magnitude below it. In the future this difference 
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Figure 10: Scatter plots of M^. — M^o versus (av) for different numbers of 
coannihilating N: 1 (orange diamonds), 2 (blue circles) and 3 (red squares). 



will remain. The smallest ctsi that could be excluded by XENONIT is about 
10~^° pb for coannihilating neutrinos but 3 x 10"^^ pb for 3 coannihilating 
neutrinos. Coannihilations, therefore, may have an important impact in the 
direct detection prospects of dark matter within this model. 

What about indirect detection? Well, it turns out that Ni coannihilations 
may have even more significant implications for the indirect detection of dark 
matter, as they allow to reconcile a thermal freeze-out with large values of {(Jv), 
and therefore with enhanced indirect detection signals. Figure [TUl shows (av) as 
a function of M^. — M^o for different number of coannihilating neutrinos. We 
see that the higher the number of coannihilating neutrinos the larger {av) can 
be. In fact, {av) goes from a maximum of 2 x 10~^^cm'^s~^ for 1 degenerate 
neutrino to a maximum of 10~^*cm'^s~^ for 3 degenerate neutrinos. It is also 
clear that {av) decreases as the N.i-H'^ increases, so the largest {av) correspond 
to the smallest Ma?. — Muo . 

A comparison of the viable regions in the plane {Mho ,{av)) is presented in 
figure [TTl The yellow region (the most internal one) shows the viable region in 
the absence of coannihilations with the right-handed neutrinos. That is, it is 
the result for the inert doublet model. Notice that in that case the minimum 
viable mass is about M^a ^ 500 GeV and that {av) lies slightly above its 
thermal value, approximately between 4 x 10~^^cm'^s^^ and 8 x 10~^^cm'^s~^. 
As we include coannihilating neutrinos, the minimum value of Mj^o decreases 
and the range of variation of {av) increases. For 1 coannihilating neutrino, 
the viable region starts at Af^o ^ 300 GeV and the maximum value of {av) 
reaches almost 3 x 10~^^cm^s~^. If all three neutrinos coannihilate with 
one can obtain viable models already at M^o ~ 100 GeV and {av) might be 
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Figure 11: Regions in the plance {Mjjo ,{av)) that are consistent with the dark 
matter constraint for different numbers of coannihilating N. 



above 10~^'*cm^s~"'^. More generally, this figure demonstrates that there is no 
incompatibility between a large value of (av) (say of order few x 10~^^cm^s^^) 
and the idea of a thermal freeze-out. Both can be reconciled via coannihilations 
in the early Universe. 

Figure [TT] also shows the current bound on (crw) from Fermi-LAT [3T] . It 
excludes the low mass region, M^o < 250 GeV, and also the models with the 
largest values of (av) over the entire mass range we consider. There are still, 
however, sizable regions where the coannihilation effects are important that are 
not constrained by present data. Future indirect detection data will certainly 
probe this region further, either improving the constraints or finding evidence 
of its existence. 

Summarizing, coannihilations with right-handed neutrinos may play a very 
important role in the dark matter phenomenology of the radiative seesaw model. 
They open up new viable regions, change the expected mass spectrum, and 
modify the prospects for the direct and the indirect detection of dark matter. 



6 Discussion 

We have throughout this paper assumed a specific structure for the matrix of 
neutrino Yukawa couplings, equation 1^. One may wonder, therefore, how our 
results depend on that particular choice. From the discussion on coannihilations, 
section [31 we know that for Ni coannihilations to increase the relic density it 
must happen that (Thoho ^ o'HONi, '^NiNi- So, all we need to reproduce our 
results is to ensure that hai are such that the above conditions on the cross 
sections are satisfied, and it turns out that they usually are. In fact, we have 
repeated the scan over the parameter space of the model (see section [5]) but 



15 



allowing the 9 parameters hai to vary randomly without being subject to any 
constraint related to neutrino masses. The results so obtained are indistinguish- 
able from those presented here. We can thus conclude that the effects discussed 
in this paper do not depend on the choice made in equation © and, more gen- 
erally, that they are not that sensitive to the experimental data on neutrino 
masses. 

The discovery potential of the inert doublet at the LHC have been ex- 
plored in the dilepton [351 [33], trilepton [5^, and multilepton [3S] channels 
plus missing energy. All these studies have focused on the light mass regime 
Mho < 80 GeV [SB]. The high mass regime. Mho > 500 GeV, features small 
mass splittings between the odd scalars[T8l [28] and is out of the LHC reach 
because of the small cross sections and rather soft leptons and jets. From fig- 
ures [7] and [5] we can see that now we can have substantially lower masses and 
much larger splittings. However, due to their large backgrounds, we still expect 
the signals for Mho > 200 GeV to be difficult to detect. The main effect of 
the coannihilating neutrinos at colliders in this scenario with scalar dark matter 
could be the modification of the cascades leading to H^. If some of the Yukawa 
couplings hai are competitive with the gauge couplings, and the spectrum is 
such that Mh± > M^. {Mj^o > MjvJ, then the leptonic signals (missing en- 
ergy) could be increased through processes like — >■ l^vH^ — > vi/H^) 
mediated by iV^. A more detailed study of these signals and of their possible 
relevance for LHC searches is left for future work. 

As we have seen, in the radiative seesaw model, Ni coannihilations allow to 
boost (av) well above the thermal value of 3 x 10~^®cm'^s^^, enhancing in a 
significant way the indirect detection signals in this model. The idea of using 
coannihilations to increase the indirect detection rate is certainly not new; it was 
examined, for example, in [26]. But it seems to us that, in spite of its possible 
relevance to current indirect detection data, it has not received much attention 
lately. For that reason, we would like to emphasize that this setup is quite ver- 
satile and can be applied within different models of dark matter. The necessary 
ingredients are just two: a dark matter particle with a large self-annihilation 
cross section and a number of coannihilating partners with weaker interactions. 
The required number of coannihilating particles and their mass splitting is cer- 
tainly model dependent but it can easily be calculated on a case by case basis. 
Notice that, among the different possibilities that have been considered to in- 
crease the dark matter indirect detection rate (gravitational clustering of dark 
matter, non-thermal production of dark matter, non-standard freeze-out, etc.), 
coannihilations seems to be the simplest one, allowing to preserve the idea of a 
thermal freeze-out and making definite predictions about the particle spectrum 
of physics beyond the Standard Model. 

7 Conclusions 

The radiative seesaw model is a simple extension of the SM that can account 
for neutrino masses and explain the dark matter. Within this model, we have 
considered the case in which the dark matter candidate is a scalar, and have 
studied the possible effect on the dark matter density of coannihilations with the 
right-handed neutrinos Ni. The dependence of these coannihilation effects with 
the relevant parameters of the model -the mass splitting, the neutrino Yukawa 
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couplings, and the number of coannihilating neutrinos was analyzed in detail. 
The net effect of the coannihilations with the right-handed neutrinos has been 
found to be an increase in the relic density, allowing to satisfy the dark mat- 
ter constraint for smaller values of M^o. Its lower bound, in fact, moves from 
about 500 GeV (without Ni coannihilations) to about 100 GeV (with all three 
Ni coannihilating simultaneously). Consequently, the dark matter phenomenol- 
ogy of this model can be substantially modified and may differ in important 
ways from that of the inert doublet model. Making use of a random scan over 
the entire parameter space of this model, we have investigated the new viable 
regions that arc obtained when coannihilations arc taken into account. Such 
regions were shown to feature bigger values for the mass splitting between the 
scalar particles and improved prospects for the direct and the indirect detec- 
tion of dark matter. Remarkably, one can obtain values of {av) as large as 
10~^^cm^s~^ (well above that normally associated with a thermal relic), with 
important implications for the indirect detection of dark matter. Finally, we 
emphasized that one could use coannihilations to reconcile large values of (av) 
with a thermal freeze-out also within other models of dark matter. 
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